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Abstract 

The Groverian entanglement measure, G{ip), is applied to characterize a variety of pure quantum 
states \tjj) of multiple qubits. The Groverian measure is calculated analytically for certain states of 
high symmetry, while for arbitrary states it is evaluated using a numerical procedure. In particular, 
it is calculated for the class of Greenberger-Horne-Zeilinger states, the W states as well as for 
random pure states of n qubits. The entanglement generated by Grover's algorithm is evaluated 
by calculating G{tjj) for the intermediate states that are obtained after t Grover iterations, for 
various initial states and for different sets of the marked states. 
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I. INTRODUCTION 



In the past decade it was demonstrated that for certain computational tasks, quantum 
algorithms are more efficient than classical algorithms. Two examples of this effect are pro- 
vided by Shor's factorization algorithm jl| and Grover's search algorithm 0,0]. While the 
understanding of quantum algorithms is still incomplete, there are indications that quantum 
entanglement plays a crucial role in making them powerful. Quantum algorithms generate 
entangled states that involve large numbers of qubits. To assess the role of entanglement 
in the algorithms it would be useful to develop ways to quantify it, which are based on 
operational considerations. Such considerations in the context of quantum communication 
between two parties were successfully applied to develop measures of entanglement for bi- 
partite systems 0,0,0. It is thus expected that operational considerations in the context of 
quantum computation, that involves multiple qubits, may lead to useful entanglement mea- 
sures for multi-partite systems. Recent work based on axiomatic considerations has provided 
a set of properties that entanglement measures should satisfy 

mm 

. These properties 

include the requirement that any entanglement measure should vanish for product states, 
it should be invariant under local unitary operations and should not increase as a result of 
any sequence of local operations complemented by only classical communication between the 
parties. These properties provide useful guidelines in the search for operational measures of 
entanglement based on quantum algorithms. Measures that satisfy the properties specified 
above are called entanglement monotones. 

Recently, it was shown that the entanglement of a pure quantum state \ip) of n qubits 
can be quantified by the utility of this state as an initial state for Grover's search algorithm 
Suppose that a given state is used as the initial state. However, before starting the 
search, one is allowed to perform local unitary operations on each qubit in order to maximize 
the probability, Pmax? of having a successful quantum search. It was shown that if {ip) can 
by transformed to some other state |0) by local operations and classical communication 
then Pma^i'ip) < -Pmax(0)- This observation was used in order to construct the Groverian 
entanglement measure Glip) of pure states. The Groverian measure is an entanglement 
monotone, and is equivalent to an entanglement measure proposed previously in Refs. , 
which is based on distance measures between lip) and the nearest disentangled state. 

In this paper we present an explicit representation of G (■?/') and calculate it for a variety 
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of pure entangled states using a combination of analytical and numerical methods. The 
difficulty is that the calculation of the Groverian measure is based on the maximization of 
a multi-variable function. The number of variables over which the maximization is taken 
increases with the number of qubits, ra, in the state The measure G{il)) is calculated 
analytically for various quantum states that exhibit high symmetry, such as the Greenberger- 



Horne-Zeilinger state 



12j | and the W state |13|, ll^ of n qubits. The numerical procedure is 



applied to the calculation of G{il)) for the states that are produced by Grover's algorithm 
as a function of the number of iterations, t, when the initial state is either a product state 
or an entangled state and for different sets of marked states. Finally, random states of n 
qubits are also examined. It is found that homogeneous random sampling of pure states of 
n qubits produces highly entangled states for which G{il)) approaches 1 for large n. 

The paper is organized as follows. The quantum search algorithm, using arbitrary initial 
states is described in Sec. II. The Groverian entanglement measure is presented in Sec. III. 
An explicit representation of the Groverian measure is introduced in Sec. IV. This explicit 
representation is used in Sec. V to perform analytical calculations of the Groverian measure 
for certain pure states of high symmetry. Numerical calculations of G'(^) for arbitrary pure 
states are presented in Sec. VI. The results are summarized in Sec. VII. 



II. QUANTUM SEARCH USING AN ARBITRARY INITIAL STATE 

Consider a search space D containing elements. We assume, for convenience, that 
= 2", where n is an integer. In this way, we may represent the elements of D using an 
n-qubit register containing the indices, i = 0, . . . , N — 1. We assume that a subset of r 
elements in the search space are marked, that is, they are solutions to the search problem. 
The distinction between the marked and unmarked elements can be expressed by a suitable 
function, f : D ^ {0, 1}, such that / = 1 for the marked elements, and / = for the rest. 
The search for a marked element now becomes a search for an element for which / = 1. To 
solve this problem on a classical computer one needs to evaluate / for each element, one 
by one, until a marked state is found. Thus, on average, N/2 evaluations of / are required 
and in the worst case. For a quantum computer, on which / to be evaluated coherently, 
it was shown that a sequence of unitary operations called Grover's algorithm can locate a 
marked element using only 0{\jN/r) coherent queries of / 
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To describe the operation of the quantum search algorithm we first introduce a register, 
\i) — \ii . . .in), of n qubits, and an ancilla qubit, \q), to be used in the computation. We 
also introduce a quantum oracle, a unitary operator O which functions as a black box with 
the ability to recognize solutions to the search problem. The oracle performs the following 
unitary operation on computational basis states of the register, and the ancilla, \q): 

0\i)\q)^\i)\q®f(i)), (1) 

where © denotes addition modulo 2. The oracle recognizes marked states in the sense that 
if \i) is a marked clement of the search space, namely f{i) = 1, the oracle flips the ancilla 
qubit from |0) to |1) and vice versa, while for unmarked states the ancilla is unchanged. The 
ancilla qubit is initially set to the state 

|-). = ^(|0)-|1)). (2) 
With this choice, the action of the oracle is: 

= (3) 

Thus, the only effect of the oracle is to apply a phase of —1 if \i) is a marked basis state, 
and no phase change if \i) is unmarked. Since the state of the ancilla does not change, one 
my omit it and write the action of the oracle as 0\i) — (— 

The original algorithm, as introduced by Grover, includes an initialization stage in which 
the n + 1 qubits of the register and the ancilla are prepared in the state |0)®"|0)g. Then, 
the following procedure is performed: a Hadamard gate H = -^{ l }i ) is applied on each 
qubit in the register, and the gate HX on the ancilla, where X = ( ? J ) is the NOT gate. 
The matrices are expressed in the computational basis (|0), |1)). The resulting state is: 

1^)1-),, (4) 

where 

The state r) is considered as the intial state of the algorithm. 
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Here we consider a generalized algorithm in which an arbitrary, possibly entangled state 

Af-l 
1=0 

is used as the initial state instead of |?7). The ancilla is still prepared as before, namely its 
state is \—)q. The algorithm itself consists of r repetitions of the following Grover iteration: 

1. Apply the oracle, which has the effect of rotating the marked states by a phase of tt 
radians. Since the ancilla is always in the state \—)q the effect of this operation may 
be described by the following unitary operator 

^; = E(-ip^K)(^i, (7) 

i=0 

acting only on the register. 

2. Carry out the following steps: (i) apply the Hadamard gate to each qubit in the 
register; (ii) Rotate the |00 . . . 0) state of the register by a phase of tt radians. This 
rotation is similar to step 1, except for the fact that here it is performed on a known 
state. It takes the form 

N-l 

7o^ = -2|0)(0| + 5]|i)(i|, (8) 

i=0 

where the second term on the right hand side is the identity operator, denoted by I. 
(iii) Apply the Hadamard gate again on each qubit in the register. 

The resulting operation is 

_^®n^^®n ^ _7 + 2i7®"|0)(0|i/^" = -I + 2\r]){r]\. (9) 
When this operator is apphed on the state l^') it results in the state 

N-l 

_H^nj.^m\^^ ^ J2{2a - a,)\i), (10) 

i=0 
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where 



Af-l 



1=0 



(11) 



Thus, each amphtude is rotated by vr around the average of all amplitides of the 
quantum state. 

The combined operation on the register in one Grover iteration is described by 



G 



_ jjiSin jn jj^n JIT 



(12) 



After the completion of r Grover iterations, the register is measured in the computational 



basis. The optimal number of iterations is 



r 



or, approximately for r <^ N 



N -r 



/cos-^(l - 2r/N) 



(13) 



r 



(14) 



where \_x\ is the largest integer which is smaller than x. Using the original initialization 
process, at this optimal time, a marked state can be found with almfost certainty, or more 
precisely with probability 



= 1-0 



1 



(15) 

in the sense that 



With this performance, Grover's algorithm was found to be optimal 

it is as efficient as theoretically possible |l7|. A variety of applications were developed, in 
which the algorithm is used in the solution of other problems 

00,0, HQ- The 

algorithm was also generalized by allowing an arbitrary (but fixed) unitary transformation 
to take place of the Hadamard transform and an arbitrary phase rotation instead of the tt 



mversion 



ion ^ 



24 



When a general pure state is taken as the initial state for Grover's algorithm the 
success probability is reduced. In this case, the probability P(t) to find a marked state if 
a measurement is taken after t iterations depends not only on the initial state and the 
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number of marked states, but also on the specific identity of the set of marked states. For a 
given choice of the set of maked states, the states \ip{t)) obtained after t Grover iterations, 
starting with an arbitrary initial state \iIj{0)), were calculated using recursion equations 
The optimal time to measure as well as the maximal probability of success were found to 
depend both on the initial state lip) and on the specific choice of the set of marked states. 
Of course, in a real search, the set of marked states is unknown, although it is assumed 
that the number of marked states is known. To evaluate the success probability and the 
optimal time to measure one needs to perform an average over all possible choices of the 
set of marked states. Under these conditions, the optimal time to measure was found to be 
equal to r, namely, after the same number of iterations as for the case in which the initial 

success probability Ps{ip) is reduced. It can be expressed in 
terms of the amplitudes of the initial state according to 



Thus, the success probability of the search depends only on the first moment of the distri- 
bution of the amplitudes and not on higher moments. Moreover, it does not depend on the 
number of marked states, up to a correction of order r/N j^]. 

III. THE GROVERIAN ENTANGLEMENT MEASURE 

Consider n parties sharing a pure quantum state of qubits, where each party is in 
possession of one qubit. The parties use those particular n qubits to perform a quantum 
search in the space of = 2" elements. Prior to the search, each party may perform local 
unitary operations on the qubit in their possession. After they complete the local processing 
of their qubits, all parties send (or teleport) their qubits to the search processing unit. 
The only processing available in this unit is Grover's search iterations and the subsequent 
measurement. Thus, the only way the qubits are allowed to interact is through Grover 
iterations. 

The local pre-processing can be expressed by 




(16) 



V = Ui®U2®---®Ur 



(17) 
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where Uk is an arbitrary local unitary gate acting on the kth. qubit. The initial state 
inserted into the Grover iterations is then V\il)). These local operations are chosen such 
that the success probability Pg of the algorithm will be maximized. The maximal success 
probability under these conditions will be 

Pmax(^)= max Ps(t/l ® f/n|V^)). (18) 

U\,...,Un 

It turns out that Pmax(V^) can be used to quantify the entanglement present in the state 
To make this assertion more precise, let us write Pmaxiip) in terms of the operator Uq 
representing r Grover iterations. For simplicity consider the case in which there is only a 
single marked state \m). Since the marked state is unknown, the probability Pmax(V') should 
be averaged over all possible identities of m, namely m = 0,l,...,A^ — 1. Performing this 
average we obtain that 

N-l 

P^ax(V')= max y^\{m\U^a{Ui®U2®---®Un)m\ (19) 

Ul,...,Un iV •^"^ 

m=0 

where the maximization is over all local unitary operations Ui,...,Un on the respective 
qubits. In order to evaluate Pmax(V') we recall that 

U^^\r^) = \m) + o(^], (20) 



where the second term is a small correction. Multiplying this equation by {UqY and taking 
the Hermitian conjugate gives 

{m\U-^ = {V\ + O (^-^y (21) 
Substituting into Eq. (fT^ gives, for a general state \ip), 

1^-1 / 1 \ 

P^ax(V') = ,max - y^K^lUi ® f/2 ® • ■ ■ ® UMf + 0(^]. (22) 

Uu...,Un iv WN / 

However, |?7) is a product state, so that Ul Ul ■ ■ ■ Ul\r]) is another product state. 
Therefore, the maximization in Eq. may, equivalently, be expressed by 

Pmax(V') 

max |(ei, . . . , e„|?/')| +0[—=\, (23) 

|ei,...,e„> VVA'/ 
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where the maximization now runs over all product states, |ei, . . . , e„) = |ei) ® • ■ ■ ® |e„), of 
the n qubits. 

In Ref. 01 it was shown that the maximum success probability, -Pmax(V'); can be used 
to define an entanglement measure, the Groverian entanglement, for arbitrary pure multiple 
qubit states. The Groverian entanglement of a state \ip) is given by: 



G(^) = Vl-i'maxW. (24) 

Since Pmax('?/') takes values in the range < Pmax(^) < 1, it follows that < GiTp) < 1- 
It is clear from the definition that all the states \ip) that can be reached from \ip) by local 
unitary operations share the same measure, given by G{tp). It is also easy to see that for 



all product states G(tti) = 0. In Ref. Ill it was shown that the Groverian entanglement 

uu 

measure is closely related to an entanglement measure introduced previously in Refs. j7|, |9|| . 
This relation was used in order to demonstrate that G{iIj) is an entanglement monotone, 
namely it cannot be increased by local operations and classical communication. Therefore, 
the Groverian measure is a good entanglement measure for pure quantum states of multiple 
qubits. 

IV. EXPLICIT REPRESENTATION OF THE GROVERIAN MEASURE 

Consider a pure state 

N-l 



J^fliK) (25) 



i=0 



of n qubits, where 



\i) = |zi,Z2, ■■■,in) (26) 

and ik is the kth most significant bit of the binary integer i. The Groverian entanglement 
measure G^ip) is given by Eq. (jH)), where Pmax(V') is given by Eq. (j^ . To obtain an 
explicit formula for G{ip) we will consider a general product state of n qubits 

|e) = |ei)(g)---® |e„). (27) 
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The single qubit states can be represented by 



\ek) = cos9k\0)k + e"^'' smek\l)k, (28) 

where k — 1,2, ... ,n, and global phases are ignored. Note that our angle < < 7r/2 
is a half of the angle 9 used in the Bloch sphere representation, while (pk is in the range 
0<ipk< 2n. 

The product state takes the form 



|ei) ® • • • ® |en) = cos6'i . . .cos6'„|0. . .0) + (29) 

+ cos^i...e^^"sin^„|0...01) 

+ •••+ (30) 

+ e^'^isin^i...e^'^"sin^„|l...l). 



Therefore, the product state |e) can by written in the form 



N-l 



^) = J2c^\^), (31) 



i=0 



where the coefficient of the basis state \i) = 12, ... , in) is 



n 



Q = J](cos^fe)^'=(e*'^'=sin^fc)*^ (32) 
fe=i 

where — 1 — ik- The overlap between the given state and a product state |e) is defined 

by 

N-l 

(elV') = ^ a, c*, (33) 

where c* is the complex conjugate of q. The calculation of the Groverian measure involves 
the the maximization of the function 

P{9i,...,9n,^pi,...,^Pn,tp) = |(ei,...,e„|V')|^ (34) 

with respect to the variables 9k, fk, k = l,...,n. Note that for 9k = rc/A and ipk = 0, 
k — l,...,n, this function coincides with the success probability of Grover's algorithm 
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starting with the initial state IV'), given by Psi'ip) = The maximal probability of 

success can now be written as 

PmM) = max P(6'i,...,6'„,v?i,...,V5n, (35) 

9l,...,8„,ipi,...,ipn 

up to a correction term of order l/-\/iV, and the maximization is taken in the range < 
6k < tt/2 and < v^fc < 2n. 

The case of states {ip) in which all the amplitudes a, are real is simpler. For such states, the 
product state |ei, . . . , e„) for which the maximum in Eq. (jH^j) is obtained has real amplitueds 
as well, namely all the angles = or vr. Therefore, in this case the maximization over (pk 
is reduced to a discrete maximization over the binary choice of exp{i(fk) = ±1. This term 
can be removed by doubling the range of 6k to —7i/2 < 6k < 71/2, thus allowing sin^^ to be 
both positive and negative, for the same value of cos 6k. Thus, for states \ip) in which all 
the amplitudes aj are real 

-Pmax(V^) = max P(6'i, . . . ,6'n,V'), (36) 

6l,...,6n 

where the maximization is over the range — 7r/2 < 6k < /2 and ifk = 0,k = l,...,n. 

V. ANALYTICAL CALCULATIONS OF THE GROVERIAN MEASURE 

Having found an explicit expression for the Groverian entanglement measure, we can 
now use it in order to characterize certain quantum states that are encountered in various 
contexts of quantum computation and communication. 

A. Two-qubit states 

Consider a state {ip) of two qubits in which all the amplitudes are real. According to 
Eq. (jHBj) we can write the overlap of \ip) with a tensor product of two single qubit states as 

P{6i, 62, ip) = (aoo cos 61 cos 62 + aoi cos 61 sin 62 + aio sin 61 cos 62 + an sin 61 sin 62Y. (37) 
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Using standard trigonometric identities we obtain 



r>ia a i\ /'^oo — ctoo + «ii ^ Oio + ctoi . ^ , ctio — ^oi . „ \ 
"(Pp, fc'm, V^) = I ^ COS fc'p H ^ COS Orn H ^ sm Op H ^ sm Orn I , 

(38) 

where 



Om = 01 - 02. (39) 

To obtain Pmax(V') one needs to maximize P{0i,02,i') with respect to 0i and 02, or equiva- 
lently, with respect to 0p and ft^. This maximization is done by solving 



dP{0p,0^,^) 



d0p ° 



dP{0p,0m,^P) 
d0m 

The maximun is found for 6p and 0^ that satisfy 



aoo — flu 

cost^p — 



0. (40) 



a/ {aw + aoiY + (aoo - «ii)^ 

aoo + ctii /--.s 
cos fc'm = =■ (41) 

V (aio - aoi)2 + (aoo + anV 



Inserting these values into Eq. (|HHj) gives rise to 



i'max(^) = ^ (l + Vl-4|detD|2) , (42) 
where the matrix D takes is given by 



D = I 1 . (43) 

aio ail 

Eq. ()42j) is valid also in case that the amplitudes are complex, however in this case the direct 
maximization is tedius. This result can be shown using the Schmidt decomposition and 
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the fact that PmaxlV^) is equal to the square of the maximal Schmidt coefficient [30|, l31 |. 



Therefore, the von Neumann entropy, 5*, of the reduced density matrix can be expressed by 



S = -Pmax ln2 Pmax " (1 - ^max) ln2(l " Pmax)- (44) 

Consider a generalized Bell state of the form 

1^) = aoo|00 > +aii|ll > . (45) 
Inserting these amplitudes into Eq. (jl^ we obtain that 



Pmax(^) = max(|aooP, laiiH. (46) 

Therefore, the Bell states |0±), for which |aoo| = |oii| = as well as the two other Bell 

states, are characterized by Pmax = 

In case that I?/)) is a product state, its amplitudes can be expressed by aij = biCj, i, j = 0, 1, 
where bi and cj are the amplitudes of the two single qubit states that form the state {ip)- 
Plugging this product into Eq. it is easy to see that for product states Pmax = 1- 

B. Multiple-qubit states 

1. Generalized GHZ states: 

The GHZ state of n qubits takes the form 

I^GHz) = ^(|0...0) + 11...1)). (47) 

This is a generalization of the Bell state 10+) to systems of more than two entangled qubits. 
These states can be further generalized to a continuous class of states of the form 

1^) =ao|0...0) +a;v-i|l...l), (48) 
where |aoP + |aAr_ip = 1. In this case the overlap function takes the form 

(n n \ 2 

k=l k=l / 
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To obtain Pmaxl"^) we solve the equations dP{9i, . . . ,9n,4')/d9k = and 
dP{9i, . . . ,9n, ip)/dipk = for A; = 1, . . . , ra. The solution is 9k = 0(7r/2), A; = 1, . . . , n, when 
|ao| is larger (smaller) than |aAr_i|. Therefore, Pmax(V^) = niax{|aoP, |a7v-iP} and G{ip) can 
be obtained from Eq. (j24j) . Its maximal value is obtained when |ao| = |aAr-i| = l/"\/2, 
where G(^/'ghz) = l/v^5 namely independent of the number of qubits. 

2. The W states: 

The W state of n qubits is the symmetrical superpositions of all the basis states which 
include a single qubit in the 1 state and all the other qubits in the zero state. It takes the 
form 

n 

l^w) = 7=E|2'"')- (50) 

^ k=i 

In the case of two qubits it coincides with the Bell state |^/'+) = "^(|01) + |10)), while in 
the case of three qubits it takes the form -ig(|001) + |010) + |100)). The overlap function for 
these states take the form 

2 



P{9^,...,9n) = ^iY,sin9k\[cos9k\ • (51) 

\A:=1 k'^k J 



'^k 

Taking derivatives with respect to the ^^'s we obtain that the maximal value is obtained at 
sin 9k = '^1 \/n, cos 9k = y^l — 1/n, k = 1, . . . , n. The maximal value is 



1 ^ 



n 



Pmax(^w) = 1 - - , (52) 



which converges to 1/e as the number of qubits increases. The Groverian measure then 
converges to G^ip^^) = a/1 — 1/e at n ^ oo, which is higher than the value for the GHZ 
states but lower than 1. We will now examine a class of strongly entangled states, in which 
the Groverian measure converges to 1 at n — > cxo. 

3. The balanced states 

Consider a state of an even number, n, of qubits, that consists of an equal superposition 
of all the balanced basis states, namely of all those basis states in which the number of O's 
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is equal to the number of I's. This state takes the form 



^ ( 1 0... 01... 1) + Permutations). (53) 



K 



where the binomial coefficient 



K=( ^ I (54) 
n/2 

is equal to the number of different permutations of n bits, where n/2 of them are in the 

state and the other n/2 are in the 1 state. The overlap function for this state is 



P{9i, . . . ,9n,ip) = Vk{cos 9i - ■ • cos 9n/2 sin ■ ■ ■ sin 6'„ + Permutations). (55) 

Taking the derivatives and solving for dP{9i, . . . , '^)/d9}. — 0, we obtain that the maxi- 
mum is at = 7r/4, A; = 1, . . . , n, and 

Pmax(V^) = K (56) 

In the limit of n — > 00 one can use the Stirling approximation, that gives rise to 

Pma.(^)^\/-^. (57) 

Since Pmax('0) converges to zero, the Groverian measure G('0) of the balanced states con- 
verges to 1 as n — > 00. 

VI. NUMERICAL CALCULATIONS OF THE GROVERIAN MEASURE 
A. The numerical procedure 

In order to cvahiatc the Groverian entanglement measure of a given pure state of n 
qubits one has to find the maximal overlap of this state with any product state with the same 
number of qubits. The overlap is given by the function P(6'i, . . . ,9n, (pi, ■ ■ ■ , fn, V')- This is 
thus a maximization problem in a 2n-dimensional space (which is reduced to n dimensions 
if all the amplitudes are real). In general, this problem cannot be solved analytically, and 
therefore the use of numerical calculations is essential. 
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The vector space in which the maximization is performed is 



f={9i,...,9n,(pi,...,(Pn) (58) 

for the general case in which the amphtudes are complex [in the special case of real ampli- 
tudes it is reduced to r = (^i, . . . , ^„)]. The optimization is done using the steepest descend 
method, namely 

dv -* 

-^cvp{e,,...,en,ip,,...,iPn,ij), (59) 

where C > is a constant and V = {d/d9i, . . . , d/dOn, d/d(pi, . . . , d/d(pn)- 
For a given initial point, r, the steepest descent method converge to a nearby local maximum. 
In order to obtain the global maximum, the calculation is repeated with a large number of 
random initial points. The largest value among all the local maxima that are reached is 
then picked as the numerical result for Pmax- 

In general, the function P{9i, . . . ,9n,(fi, . . . ,(fn,ip) consists of a sum of products of 
trigonometric functions. Each product includes n sine or cosine functions and is thus pe- 
riodic in all directions. The number of terms in the sum increases exponentially with n, 
making the maximization problem more difficult. We find that the number of local minima 
also increases with n. Thus, as n increases one needs more calls to the steepest descent 
program with random initial points. 

To exemplify the effect of the maximization process we consider the family of states given 

by 

= ar^\rj) + oghz | V'ghz) , (60) 

where < cghz < 1 and since the two states are not orthogonal we obtain from the 
normalization condition that 

a, = -y/2/N oghz + -^1 - (1 - 2/N)al^^. (61) 

In Fig. 1 we present the success probability Ps{ip) (dashed line) of Grover's search with the 
initial state as well as Pmax('0) (solid line) obtained from the numerical procedure, as a 
function of Oghz) 12 qubits. The two functions follow the same path, decreasing linearly 
as Oq^i increases, up to Oghz — 0.65. This means that for the states in this range the 
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success probability cannot be increased by local unitary operations. As the two functions 
depart, Pmax(^) starts to increase while Ps{ip) continues to decrease. The size of the gap 
between the two functions represents the effect of the maximization process. 

In Fig. 2 we show the Groverian measure of the generalized GHZ states [Eq. (j48|) ] of 12 
qubits as a function of |aop. The numerical results (o) coincide with the analytical results 
(solid line). The graph is symmetric around |aop = 1/2, where the largest value of G{ip) is 
obtained. At this point, for real and positive amplitudes the state {ipcRz) is obtained. 



B. Entanglement during the operation of Grover's algorithm: 

There are indications that entanglement plays an important role in making quantum 
algorithms more efficient than their classical counterparts. Therefore, it is interesting to 
see how entanglement is generated during the operation of quantum algorithms. Grover's 
algorithm is particularly suitable for this study because it consists of a large number of 
iterations of the same set of operations. Furthermore, for any given initial state \tp{0)), 
the amplitudes of the state \ip(t)), obtained after t Grover iterations, can be calculated 



analytically 
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In Fig. 3 we present the Groverian entanglement measure of the states \ip{t)), obtained 
after t = 0, 1, . . . , 50 iterations of Grover's algorithm with n = 12 qubits for the initial state 
is \ip{0)) = \ri) (dashed line). It increases with the time until it reaches its highest value 
of after 25 iterations. Then it follows the same path downwards reaching zero value 

after 50 iterations, where the marked state is reached. The case in which the initial state is 
entangled is also considered. To this end we construct the state 



l^even) + OoddlV'odd) (62) 

where \ipeven) is a normalized superposition of all the A^/2 basis states that include an 
even number of I's, while iV^odd) is a normalized superposition of all the other basis states, 
that include an odd number of I's. Apart from the initial state \ri) that is obtained for 
Oeven = l/"\/2, wc also show the results for Oeven = 0.984 (dotted line), 0.994 (dashed-dotted 
line) and 1 (solid line). The last state in this list is the one in which the amplitudes of all 
the basis states that have an odd number of I's vanish. This state can be obtained from the 
GHZ state by applying the Hadamard transform on all the qubits. 



17 



In Fig. 4 we show G(V') for the states \i>{t)) obtained after t Grover iterations where 
|'0(O)) = \ri) for 12 qubits and two marked states. In case that the marked states are |0) 
and — 1) (sohd hne) the resulting state after r iterations is IV'ghz)- In case that the two 
marked states are |0) and |1) (dashed hne) the resulting state is a superposition of these 
two states, which is not entangled. Interestingly, during the first t/2 iterations these two 
functions nearly coincide. 

In Fig. 5 we present the Groverian measure of the states obtained after t Grover 

iterations, where |'0(O)) = \r}) for 12 qubits and 12 marked states. When the set of marked 
states consists of |00 . . . 01), |00 . . . 10), . . . , |10 . . . 0), namely the basis states of non-zero 
amplitudes in the W state (solid line), the register approaches the W state after 14 iterations. 
When the marked states are i = 0, . . . ,11 (dashed line) the resulting state is entangled, 
but exhibits a smaller value of G{ip). 

C. Random states of n qubits 

How entangled is a randomly chosen pure quantum state of n qubits? To study this 
question we pick random pure states and evaluate their Groverian entanglement measure 

using the numerical procedure. 
Consider a random pure state 



of n qubits, where = |ai|e'"\ Such random states can be obtained as follows. In the 
first step one draws N numbers, Oj, i — 0, . . . , N — 1, independently from a Gaussian 
distribution centered at 0, with a standard deviation a — 1. These numbers are then 
normalized according to 



The moduli |aj|, i = 0, . . . , N — 1 are then obtained as the absolute values of the a^'s. 
The arguments ctj, i — 0, . . . , N — 1 are drawn from a homogeneous distribution in the 
range [0, 27r) . Due to the properties of the Gaussian distribution, the resulting states are 



N-l 



(63) 



i=0 




(64) 
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distributed randomly and isotropically in the Hilbert space that consists of all the pure 
states of n qubits. 

From the normalization condition we obtain that the second moment of the distribution 
of the moduli of the amplitudes, 

H' = iv5Z'^^''' ^^^^ 

satisfies |ap = namely, its square root is of order l/\fN. In the case of the randomly 
chosen states, due to the random phases, the average amplitude a [given by Eq. (fTT|) ] will be 
of the order of a = Therefore, the probability of success -Ps(V') of Grover's algorithm 

using a random state as the initial state will vanish like Our analysis of these states 

using the numerical procedure shows that for the random states the effect of the optimization 
procedure is negligible. As a result, Pmax(V^) also vanishes like Thus, the entanglement 
measure G{il)) — 1 as the number of qubits increases indicating that the vast majority of 
the states of multiple qubits are very highly entangled. 



VII. SUMMARY AND DISCUSSION 



The Groverian entanglement measure was applied to characterize a variety of pure quan- 
tum states of multiple qubits. For certain classes of states of high symmetry, the Groverian 
measure was calculated analytically. In order to evaluate it for arbitrary states with complex 
amplitudes, a numerical minimization procedure, based on the steepest descent algorithm 
was developed. It was used in order to evaluate the amount of entanglement generated by 
Grover's algorithm for different initial states and for different sets of marked states. It was 
also shown that the typical pure states of n qubits obtained by random sampling are highly 
entangled, namely for these states G{il)) = 1 up to corrections of order 

In recent years several entanglement monotones were proposed as measures of entangle- 
ment of multiple qubits 3, 0, 1^ 1^ 1^, 1^ • Unlike the case of pure states of two qubits, 
in which the von Neumann entropy provides a complete characterization of the entangle- 
ment, multiple qubit states support a large number of different measures. As a result there 
is no meaningful way to compare between the different measures. It seems that the issue 
of what measure is relevant depends on the specific physical or operational context. The 
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actual evaluation of entanglement measures turns out to be a difficult computational prob- 
lem. This is due to the fact that these measures are typically defined as an extremum of 
some multi- variable function. A singular result in this context is the explicit formula for the 
entanglement of formation of mixed states of two qubits, obtained in Refs. 

The Groverian measure was originally introduced for the case of a single marked state. 
Recently, it was shown that the same result for Pmax(V') is obtained in the case of r marked 
states, up to a correction of order r /N [28], thus removing the restriction of a single marked 
state. Grover's algorithm can be generalized by re jla cing the Hadamard transform by an 



arbitrary unitary operator U on the n qubits |2a,i2j,i2a,i2^ [ . Using the generalized algorithm 
to evaluate -Pmax(V'); one observes that the same result is obtained as long as the operator U 
is a tensor product of n unitary single qubit operators. In case that the operator U creates 
entanglement between the qubits in the register it cannot be used to evaluate the Groverian 
measure. 

Interestingly, the Groverian measure coincides for pure states with the measure proposed 
by Vedral, Plenio, Rippin and Knight 0, 0]- While that measure also applies in the case of 
mixed states of multiple qubits, we have not been able to extend the operational interpreta- 
tion of the Groverian measure beyond the case of pure states [11]. It would be interesting to 
examine the relevance of the Groverian measure to other quantum algorithms such as Shor's 
algorithm. Further studies of the Groverian measure and related concepts will hopefully 
contribute to the understanding of the role of entanglement in making quantum algorithms 
powerful. 
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FIG. 1: The success probability Ps{ip) of Grover's algorithm using {ip) as the initial state (dashed 
line), and the maximal success probability Pmax(V') (solid line), for = a^\r]) + acHzlV'GHz)) 
as a function of Oqjjz- The two functions coincide up to Oq^iz — 0.65, while above this point a 
gap appears. The gap broadens as the GHZ state is approached, demonstrating the effect of the 
maximization by the local pre-processing. 

FIG. 2: Analytical results (solid line) and numerical results (o) for the Groverian entanglement 
measure of generalized GHZ states [Eq. H48() ] as a function of |aop. The largest value of G(V') = 
l/\/2 is obtained for |aoP = 1/2, where the state coincides with the GHZ state. 

FIG. 3: The Groverian entanglement measure G{'4){t)) for states that are generated by Grover's 
algorithm as a function of the number of iterations t for n = 12 and one marked state, where 
the initial state is |^(0)). The curves were obtained for the states |^(0)) given by Eq. with 
«even = l/v^, namely the |t/) state (dashed line), 0.984 (dotted line), 0.994 (dashed-dotted line) 
and 1, namely the state -ff'^"|'0GHz) (solid line). 

FIG. 4: The Groverian measure for the states |V'(^)) obtained after t Grover iterations where 
1-0(0)) = |r/), for 12 qubits and two marked states. In case that |0) and — 1) are marked (solid 
line) the resulting state after r iterations is IV'ghz)- If, instead, the marked states are |0) and |1) 
(dashed line), the register approaches a superposition of these two states, which is not entangled. 

FIG. 5: The measure G{'ij){t)) for states that are generated by Grover's algorithm as a function of 
the number of iterations t for n = 12 and 12 marked states, where the initial state is |r/). When 
the marked states are |00 . . . 01), |00 . . . 10), . . . , |10 . . . 0) (solid line), the state obtained after 14 
iterations is the W state (up to a tiny correction due to the discrete nature of the iterations). 
When the marked states are |i), i = 0, . . . , 11 (dashed line), the resulting state is not as strongly 
entangled. 
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